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Abstract 

We study the Cauchy problem for the Whitham modulation equations for monotone increasing 
smooth initial data. The Whitham equations are a collection of one-dimensional quasi-linear 
hyperbolic systems. This collection of systems is enumerated by the genus g=0,l,2, ... of the 
corresponding hyperelliptic Riemann surface. Each of these systems can be integrated by the so 
called hodograph transform introduced by Tsarev. A key step in the integration process is the 
solution of the Tsarev linear overdetermined system. For each g > 0, we construct the unique 
solution of the Tsarev system, which matches the genus g+1 and g — 1 solutions on the transition 
boundaries. Next we characterize initial data such that the solution of the Whitham equations 
has genus g < N , N > 0, for all real t > and x. 



1 Introduction 

The Whitham equations are a collection of quasi- linear hyperbolic systems of the form |l| , (|] , |^ 

- Aj(ui,U2, . . . ,W2g+i)-^ = 0, x,t,Ui£]R, i ^ 1, ...,2g + 1, 5 = 0,1,2,..., (1.1) 
with the ordering iti > U2 > • • • > U2g+i. For a given g the system ( |1.1| ) is called g-phase Whitham 



equations. Because of the diagonal form of systems (1.1) the dependent variables ui > U2 > • • • > U2g+i 



are called Riemann invariants. For 5 > the speeds Xi{ui, U2, ■ ■ ■ , M2g+i), i = Ij 2, . . . , 25 + 1, depend 
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through Ml, ... , U2£(+i on complete hyperelhptic integrals of genus g. For this reason the g-phase system 
is also called genus g system. The zero-phase Whitham equation has the form 



3u _ ^^du _ ^ 
dt dx ' 



(1.2) 



where we use the notation ui = u. In this paper we study the initial value problem of the Whitham 
equations for monotone increasing smooth (C°°) initial data u{x,t = 0) = uq{x), where the range of 
^0(2;), X e ]R, is the interval (a, b) and —00 < a < b < +00. 

The initial value problem consists of the following. We consider the evolution on the x — u plane of 
the initial curve u(x,t = 0) = Uo{x) according to the zero-phase equation (1.2). The solution u{x,t) 



of (1.2) with the initial data uq{x) is given by the characteristic equation 



X = —6tu + f(u) 



(1.3) 



where f{u) is the inverse function of ^0(2;). The solution u{x,t) in (O) is globally well defined only 



for < t < tQ, where to = ^ min„gifi[/'(M)] is the time of gradient catastrophe of (1.3). Near the 
point of gradient catastrophe and for a short time i > tg, the evolving curve is given by a multivalued 
function with three branches b > ui{x,t) > U2{x,t) > U3{x,t) > a, which evolve according to the 
one-phase Whitham equations. 



u(x,t,) 



=0 




x"(t). 




g(x,t)=0 



(x,t)=0 



(X„to) 
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Figure 1.1: on picture a) the dashed line represents the formal solution of the zero-phase 
equation the continuous line represents the solution of the one-phase equations. The solution 
{ui{x, t) , U2{x, t),u^(x, t)) of the one-phase equations and the position of the boundaries x^ (t), x^{t) 
are to be determined from the condition u{x~{t),t) = ui{x~ {t),t), Ux{x~ (t),t) = uix{x^ {t),t), 
u{x~^{t),t) — U3{x^{t),t), Ux{x^ {t),t) ~ U3x{x^ {t),t), where u{x,t) is the solution of the zero-phase 
equation. Picture b) represents the functions x^{t) and x^{t) on the x ~ t plane. 
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Outside the multivalued region the solution is given by the zero-phase solution u{x,t) defined 
in ( |l.3| ). On the phase transition boundary the zero-phase solution and the one-phase solution are 
attached C-^-smoothly (see Fig |l.l| ). 

Since the Whitham equations are hyperbolic, other points of gradient catastrophe can appear in the 
branches ui{x,t) > U2{x,t) > u^^x^t) themselves or in u{x,t). 

In general, for t > to, the evolving curve is given by a multivalued function with and odd number 
of branches b > ui{x,t) > U2{x,t) > ■ ■ ■ > U2g+i{x,t) > a, g > 0. These branches evolve according 
to the g-phase Whitham equations. The 5-phase solutions for different g must be glued together 



in order to produce a C"'^-smooth curve in the {x,u) plane evolving smoothly with t (see Fig. 1.2). 
The initial value problem of the Whitham equations is to determine, for almost all t > and x, 
the phase g{x,t) > and the corresponding branches ui{x,t) > U2{x,t) > ••• > U2g+i{x,t) from 
the initial data x — f{u)\t=Q. For example for the initial data x — xq + 6tQU + {u — uq)^, to ^ 0, 
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u,(x) 


g=l 


g5_ 

u(x) 
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Figure 1.2: the plot on the x — u plane of the evolving multivalued curve for fixed i > 0. 



the solution of the Whitham equations is of genus at most equal to one [|l2| |^ . It is of genus one 
inside the cusp — 12\/3(t — tQ)i < x — xq < |y^(i — ^o)^ (l3); it is of genus zero outside this cusp. 
The point x ^ xa,t = to in the point of gradient catastrophe of the zero-phase solution. The curve 
x^{t) = xo — 12\/3(t — io)^, i > to, describes the locus of points where ui{x,t) = U2{x,t). The curve 
x^{t) = xq + f y^(^ ^ ^0)^ , t > to, describes the locus of points where U2{x, t) = U3(x, t). For generic 
initial data it is not known whether the genus of the solution of the Whitham equations is bounded. 
We say that the Cauchy problem for the Whitham equations has a global solution if the genus g < 00 
for all X and t > 0. 

Using the geometric-Hamiltonian structure Q| of the Whitham equations, Tsarev Q showed that 
these equations can be locally integrated by a generalization of the method of characteristic. Namely 
he proved that if the fimctions Wi ~ Wi{ui, U2, ■ ■ ■ , U2g+i), * — 1; • • ■ j 2g -I- 1, solve the linear over- 
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determined system 

dwi 1 d\i 



duj Xi — Xj duj 



[w,-wj], i,j = 1,2,..., 25 + 1, (1.4) 



where Xi — Xi{ui, U2, . . . , U2g+i), i = 1, . . . ,2g + 1, are the speeds in (LI), then the solution u{x, t) 



{ui{x,t),U2{x,t), . . . ,U2g+i{x,t)) of the so cahed hodograph transform 

X — —Xi{u) t + Wi{u) i = 1, . . . , 2g + 1 , (1-5) 



satisfies system (LI). Conversely, any solution {ui{x,t),U2{x,t), . . . ,U2g+i{x,t)) of (LI) can be ob- 
tained in this way. 

Tsarev theorem relies on two factors: 

a) the existence of a solution of the linear over-determined system ( |l.4| ); 

b) the existence of a real solutions ui{x^t) > U2{x,t) > ••• > U2g+i(a;,i) of the hodograph 



transform (1.5) 



In this paper we solve completely problem a), namely we build for any smooth monotone increasing 
initial data x — f{u)\t=o the solution of system ( |l.4| ) foranyy > 0. The solution (t?) , i — 1,. .. ,2g+l, 
of genus g satisfies some natural boundary conditions which guarantee its uniqueness. 
Regarding problem b) we characterize initial data such that the solution of the Whitham equations 
exists only for g < N where A'' is a positive integer. Namely we show that if the initial data satisfies 
the condition 

^-5^/(^):=/(2^+i)(h)>0, l<iVeIN (1.6) 

for all real u belonging to the domain of / except one point, then the solution of the Whitham 
equations has genus at most N for any x and t > 0. For iV = 1 this result has already been proved by 
Tian §. 

The investigation of the initial value problem of the Whitham equations was initiated by Gurevich 
and Pitaevskii Q. In the case g < 1 they solved the initial value problem of system (1.1) for step-like 
initial data and studied numerically the case of cubic initial data. 

The initial value problem of the Whitham equations was deeply studied by Tian. In he con- 
structed the general solution of the Tsarev system ( |l.4| ) for < (7 < 1 and for smooth monotone 



increasing initial data. He also proved the solvability of the hodograph transform (L5) for < g < 1. 

In ^ he obtained the solution of the Tsarev system for g > for polynomial initial data. 

The equations (1.1) were found by Whitham |^ in the single phase case g — 1 and more generally 



by Flaschka, Forest and McLaughlin [|| in the multi-phase case. The Whitham equations were also 
found in when studying the zero dispersion limit of the Korteweg de Vries equation. The hyperbolic 
nature of the equations was found by Levermore j^. 
This paper is organized as follows. 

In Sec. 2 we give some background about Abelian differentials on hyperelliptic Riemann surfaces. 
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In Sec. 3 we describe the Whitham equations. 

In Sec. 4 we build the solution of the Tsarev system (1.4) for smooth monotone increasing initial data. 
In Sec. 5 we show that under the hypothesis (^]^) the hodograph transform is solvable for g < N. 
In Sec. 6 we draw the conclusions. 



2 Riemann surfaces and Abelian differentials: notations and 
definitions 

Let 

( 2g+l 1 

Sg := = {r,fi), ^J-'^ = Y\ir - Uj)> , Ui > U2 > ■ ■ ■ > U2g+i , G IR , (2.1) 

be the hyperelliptic Riemann surface of genus g > 0. We shall use the standard representation of Sg 
as a two-sheeted covering of CJP^ with cuts along the intervals 

[u2k,U2k-l], fc = 1, . . . ,3+ 1, U2g+2 = -00. 

We choose the basis {aj, (3j}j^i of the homology group Hi{Tg) so that aj lies fully on the upper sheet 
and encircles clockwise the interval [u2j, U2j~i], j — I, . . . , g, while Pj emerges on the upper sheet on 
the cut [u2j, U2j-i]j passes anti-clockwise to the lower sheet trough the cut (— oo, U2g+i] ^md return to 
the initial point through the lower sheet. 

The one-forms that are analytic on the closed Riemann surface Sg except for a finite number of points 

are called Abelian differentials. 

We define on Sg the following differentials |l0[ : 

1) The canonical basis of holomorphic one- forms or Abelian differentials of the first kind (f>i, (f>2 ■ ■ ■ 4'g'- 



^g-l k ^ j.g-2 ^ k 

Mr)^ rV = (2-2) 



The constants are uniquely determined by the normalization conditions 

h=Sjk, j,k=l,...,g. (2.3) 

We remark that an holomorphic differential having all its a-periods equal to zero is identically zero 

2) The set cr^, fc > 0, g > 0, of Abelian differentials of the second kind with a pole of order 2fc -|- 2 at 
infinity, with asymptotic behavior 

cr^(r) = r''^^+0{r^^) dr for large r (2.4) 
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and normalized by the condition 



<-0, j = l,. 



(2.5) 



We use the notation 



CT^(r) =dpf(r), I2af{r) ^ dqS{r) g>0. 



(2.6) 



In hterature the differentials dp^{r) and dq^{r) are called quasi- momentum and quasi-energy respec- 
tively The explicit formula for the differentials u^, fc > 0, is given by the expression 



PM 



dr , P|(r) = rf+'= + a^r^+'^-i + 4rS+>'-^ ■■■ + aj^^^ , 



(2.7) 



where the coefficients — a'l{u) , u — {ui,U2, ■ ■ ■ , U2g+i), i — 1, . . . , g + k, are uniquely determined 



by (^) and (^. 

3) The Abelian differential of the third kind u)qqg{r) with first order poles at the points Q = {q,fi{q)) 
and Qo = {qo, ^J■{qo)) with residues ±1 respectively. Its periods are normalized by the relation 



(2.8) 



2.1 Riemann bilinear relations 

Let uji and UJ2 be two Abelian differentials on the Riemann surface Sg . If all the residues of lui and 
UJ2 are equal to zero, then the integrals d^^uji and d^^uj2 do not have logarithm singularities on Sg. 
If the differential loi has non zero residue, then its integral has logarithm singularities. Let s be the 
path connecting the singular points of d'^^uji. We have the following relation. 



E 



/ wi / uj2~ j ^2 1 ^1 + A(d ^uji)uj2 = 27ri Rcs[(d ^cji)w2] , (2.9) 



UJ2 / t^l 



S„-s 



where A(d"^wi) is the difference of the values of d~^uji on the two sides of the cut s and the quantity 
J2sg~s 'Rss[{d~^uJi)uj2] is the sum of the residues of the differential {d~'^LUi)uj2 on the cut surface Sg — s. 
This formula is known as the Riemann bilinear period relation p| . 
Assuming uji = ujqqg and UJ2 = (t>k h^ ( |2.9D we obtain 



2m 



,9- 



Assuming uji = ujqq^ and 0^2 — ^ppa in (2.9) we obtain 



990 ~ / '^PPo 
90 



(2.10) 



(2.11) 
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Differentiating with respect to p and q the above expression we obtain the identity 

dq[^qqo{p)\ = dp[ujpPoi<l)\ : (2-12) 

where dq and dp denote differentiation with respect to q and p respectively. From the expression ( ^.11 ) 
it foUows that ujqq„ {r) is a many- value analytic function of the variable q. The many- value character 
of LOqq^^ (r) as a function of q can be described by the equations 

f dq[LUqq,{r)] ^ Q , f (r)] = 2^Z(/.fc (r ) , fc = l,...,g, (2.13) 

In the following we mainly use the normalized differential tOzi^) which has simple poles at the 
points Q^{z) — (z,±/i(z)) with residue ±1 respectively. 

The differential i^f (f") is explicitly given by the expression 

dr J. / \ f dt fi{z) 



Lol(,r) 



/i(r) r — z 



k=l 



(2.14) 



where (fikir), fc = 1, . . . , 5, is the normalized basis of holomorphic differentials, (r) as a function of 
z, is an Abelian integral having poles of first order at the points Q^{r) ~ (r, ±/i(r)). The periods of 
this integral are obtained from the relations ( 2.13| ) 



4K(r)]=0, 



4K(r)] =47ri0fc(r), j^l,...g. 



(2.15) 



ft 



We apply the Ricmann bilinear relation (2^) to the differentials cr^{r) and iLjf{r) getting 



Res[c^f(r)d-Vf„(r)], m = 0, , 

r—oo 

4 



,9, 



2m 



- -/i(z)e™g + ^^7}T„,+i_j / 
\ j=i 



dt ii{z) 
/i(t) t — z 



(2.16) 



In the above formula e^g = 1 for m = g and zero otherwise, the coefficients have been defined in 
(2.2) and the F/'s are the coefficients of the expansion for ^ ^ 00 of 

1 



(2.17) 



We define Ffc = for fc < 0. 

Inverting (2.16) and introducing the quantities Nj{z,u) — Y^^k=i 7j —t-t ^^^"^ we obtain 



/7Vi(z,ti)\ 
^2(2, u) 

A^g(z,M) 



fi 





fo 



^ Tg Fg_l 



... o\ 

... 

fo 

fi fo/ 



4 J(3-(z) "olsJ 

4 jQ-(z) '^ly^i 



4 JQ-{z) '^g-iy^J 
_2g+l rQ+(z) g/^N ; 



(2.18) 
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where the F^'s are the coefficients of the expansion for ^ ^ oo of 



Using ( 2.18 ) the differential uj^{r) turns out to be given by the relation 

9 r9~k 



/i(r) r — z 



z, u) 



-dr . 



(2.19) 



(2.20) 



We remark that wf (r) is a multivalued function of z and it is regular at infinity. 
From the relation ( ^.18 ) we obtain the identity which will be useful later 

^i{z) = -Y^{2k-l)f,+,_J ai_^{i). 

The next proposition is also important for our subsequent considerations. 



(2.21) 



Proposition 2.1 The Abelian differentials of the second kind cr'^ir), k > 0, defined in (2.4) satisfy 
the relations 



o-^(r) = i Res wf(r)z* ^dz 



1 



-dr Res 



(2.22) 



2fc + 1 z=oo 

where wf (r) has been defined in ( 2. 14 ), t^r(^) ^•s normalized Abelian differential of the third kind 
with simple poles at the points Q^{r) — (r, ±/i(r)) with residue ±1 respectively and dr denotes differ- 
entiation with respect to r. 

Proof: the differential ReSz^oo uj^{r) z''~^ dz is normalized because uj§{r) is a normalized differ- 
ential. Furthermore 



Res 

z—oc 



wf(r)z'' ^dz = r*" ■^dr + Oir ^)dr for r ^ oo. 



Therefore Res, 



cjf(r)z -idz 



coincides with the normalized Abelian differential of the second 



kind cr^(?'). For proving the second equality in ( 2.22| ) we consider the integral in the z variable 







dAcotir) z'+-^) = I z^+-^{d^u;l{r))+i (fc - i)(c.f (r) 2^=+^) , 

Coo C*oo Coo 



where Coo is a close contour around the point at infinity. Substituting the identity dzUj^{r) = drLo9^{z) 
in the right hand side of the above relation we obtain the second relation in (2.22). □ 



3 Preliminaries on the theory of the Whitham equations 



The speeds Ai(Mi, U2, . . . , U2g+i) of the g-phase Whitham equations (1.1) are given by the ratio 

dqB{r) 



\i{u) 



dp9{r) 



, z = l,2,...,2g+l, 



(3.1) 
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where dp^{r) and dq^{r) have been defined in (E^). In the case g = 



dp°{r) - ^= , dg"(r) = dr , (3.2) 



SO that one obtains the zero-phase Whitham equation (|L 

For monotone increasing smooth initial data x — f{u)\t=o, the solution of the zero-phase equation 

(1.2) is obtained by the method of characteristic 0| and is given by the expression 

X = -6iu + /(li) . (3.3) 

The zero-phase solution is globally well defined only for < t < where to = | minug]F{,[/'(u)] is the 
time of gradient catastrophe of the solution (3.3). The breaking is cause by an inflection point in the 
initial data. For t > to we expect to have single, double and higher phase solutions. For higher genus 
the Whitham equations can be locally integrated using a generalization of the characteristic equation 

(3.3) . We have the following theorem of Tsarev g 



Theorem 3.1 // Wi{u), u — [ui, U2, . . . , M2g+i); solves the linear over-determined system 

dw. 
h 

where 

1 dX 



— ^aij{u)[wi-Wj], ij ^l,2,...,2g + l, i ^ j, (3.4) 



X < T = l,2,...,2g + l, i 7^ j, (3.5) 

\i — \j ouj 

then the solution (ui(x, t), U2{x, t), . . . , U2g+i{x, t)) of the hodograph transformation 

X = —Xi{u) t + Wi{u) i = 1, . . . ,2g + 1 , (3.6) 



satisfies system (1.1). Conversely, any solution {ui,U2, . . . ,U2g+i) of (1.1) can be obtained in this way. 

To guarantee that the 17— phase solutions for different g are attached continuously, the following 
natural boundary conditions must be imposed on Wi{ui, U2, . ■ . , U2g+i), i — 1, . . . ,2g + 1. 
When ui = ui+i, 1<1 <2g, 

wf (Wl, ...,Ul,Ui,.. . ,M2g+l) = wf^^(ui, ...,Ul,Ui,.. . ,M2g+l) (3.7) 

and for 1 < z < 2.g + 1, i ^1,1 + 1 

wf{ui, ...,Ui,Ui,.. .,U2g+l) = wf~^{ui, . . . . . . ,M2g+l). (3.8) 

The sup-scripts g and 5 — 1 in the WiS specify the corresponding genus and the hat denotes the 
variable that have been dropped. When g = 1 and U2 = we have that 

Wi(ui,U3,'U3) = /(wi) 
W2[Ui,U3,U3) = W3{ui,U3,U3) , 
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where f(u) is the initial data. Similar conditions hold true when ui = M2, namely 



W3{ui,Ui,U3) = /(U3) 
Wi{ui,Ui,U3) = W2{ui,Ui,U3). 



(3.10) 



We remark that the Xi {u) 's satisfy the boundary conditions (|3.7h3.8| ) and for g = 1 we have 

Ai(-Ui,M3,M3) = -6ui, A2(wi,U3,M3) = A3(ui,U3,U3) , 

and 

A3(ui,Ui,U3) = -6-U3, Ai(ui,Ui,U3) = A2(wi,Ul,M3). 



The solution of the boundary value problem (^^), (^-3.10) has been obtained in 
increasing analytic initial data of the form 

X = /a(w) = Co + ClU-\ h Cku'' + ... 



HM for monotone 



(3.11) 



where we assume that only a finite number of Ck is different from zero. For such initial data the Wi{uys 

oundar 



which satisfy (3.4) and the boundary conditions (3.7-3.1C) are given by the expression 

ds9(r) 



dp9{r) 



l,...,2g+l. 



The differential ds^ in ( 3.12 ) is given by 



(3.12) 



(3.13) 



and the differentials o'f.{r), k > have been defined in (p??]). 

The solution (3.6) of the 5-phase Whitham equations can also be written in an equivalent algebro- 
geometric form [pi W^ namely 



{-xdp3{r) - tdqS{r) + dsS{r)) \r= 



0, 



1,2..., 2.9 + 1 



(3.14) 



where dp^, dq^ and ds^ have been defined in ( |2.6D and ( 3.13 ) respectively. 

The solution ui > U2 > ■ ■ ■ > U2g+i of the g-phase Whitham equations (1.1) is implicitly defined as a 
function of x and t by the equations ( [3.6| ) or ( [j.l4 ) . The solution is uniquely defined only for those x 
and t such that the functions Ui{x^ t) are real and dxUi{x, t), i = 1, . . . ,2g + I, are not vanishing. 
One of the problems in the theory of the Whitham equations is to determine when (3.6) or (3.14) are 
solvable for real ui, . . . , U2g+i as functions of x and t. This problem has been solved by Tian for 5 < 1. 
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Theorem 3.2 ^ Consider a monotone increasing initial data x — f(u)\t=o. Suppose that u* is the 



inflection point of f{u) that causes the breaking of the zero-phase solution (l.i) at x — xq, t = to. 
Let be f"'{u) > in a deleted neighborhood of u = u* . Then the one-phase Whitham equations has 
a solution ui > U2 > U3 within a cusp in the x — t plane for a short time after the breaking time of 



the zero-phase solution. Furthermore this solution satisfies the boundary conditions (S.L) and (3.1C) 
on the cusp. If the initial data satisfies the condition f"'(u) > for all u except u = u* , then the 
solution of the Whitham equations exists for all t > 0. The solution is of genus one inside the cusp 
x~(i) < X < x'^{t), t > to, where x^{t) < x^{t) are two real functions satisfying the condition 
x^ita) — x^{to) — xq. The solution is of genus zero outside the cusp x~{t) < x < x~^{t), t > to. 



4 Solution of Tsarev system 



In this section we build the solution of the boundary value problem (3^), (^- 3.1C ) for monotone 
smooth initial data and we show that the solution obtained is unique. We consider initial data of the 
form X = f{u)\t=o where f{u) is a monotone increasing function. The domain of / is the interval (a, b) 
where —00 < a < b < +00, and the range of / is the real line (—00, +cx)). 



In order to obtain the solution of the boundary value problem (3^), (3.7- 3.10 ) we need the following 
technical lemma. 



Lemma 4.1 The differential ds^ir) defined in (3. IS) can be written in the form 

RB{r) 



/ 2g+l \ 

dsO{r) = 2^(r) dr-^3{r-u) + ^ 9„,*s(r;w) dr 



fiir) 



where 



^'^(r; u) 



Res 

2:— 00 



!F{z)dz 



2/z(z)(z-r) 



qk{u) 



Res 

2—00 



z3-''T{z)dz 



dr, 



k = l,...g, 



Jo 



faiO 



(4.1) 



(4.2) 



(4.3) 



2<7+l 2g+l g k 

R<^ir)=2j2duMu) n ir-ui) + Y,qk{u)J2i'^l-l)rk-iP;'_,{r), (4.4) 

k=l l = lj^k k=l 1 = 1 



and the polynomials P^{r), I > 0, have been defined in (2.1), the 's have been defined in (2.1b). 



Proof: using the second identity in ( 2.22| ), we rewrite the differential ds^{r) defined in (3.13) in the 
form 



ds^ir) = -dr f ResK(z)J^(z)]) , 

\z — oo J 



(4.5) 
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where i^riz) has been defined in (2.14) and J-{z) is the Abel transform defined in (4.3) of the analytic 
initial data (3.11). The identity (4.5) can be checked straightforward. Using the explicit expression of 
uJr(z) in ( 2.20 ) we obtain 

g k 

ds3{r) ^2dri^I{r)■9'{r;u))+Y,Qkiu)J2{2l-l)fk^laf_^{r), (4.6) 



fc=i 



1=1 



where ^^(r;u) and qkiu) have been defined in (4.2). 

2a„.*9(r;w), 2d,,^qgiu) = ^Huf,u) 



From (4.2) we get the relations 

*9(r;it) *9(ui;ii) 



r ~ Ui r ~ Ui 
and for g — Q we define ui = u and 



(4.7) 



2duqo{u) -.^^^u-u)^ fa{u). 



Using (4/7) we transform the expression for ds^lr) in (|4.6| ) to the form (4.1) 



The relation ( ^j| ) enables us to write the quantities Wi{i 



dps(r) 



I = 1, 



(3.12) in the form 



Wi{u) 



2g+l g k 

l = l,l^i k=l 1=1 



□ 

,2g + 1, in 
(4.8) 



Observe that in the formula (4.8) all the information on the initial data is contained in the functions 
qk{u). The functions qk — qk{u), k = 1, ... ,5, solve the linear over-determined system @] 



2{ui - Uj) 



d'^qkiu) dqk{u) dqk{u) 



duiduj 



qk{u,u, ...,u) 

2g+l 



dui 
(25-1)! 



, = 1, ■ ■ • ,25 + 1, 



-k 



(4.9) 



du 



where fa'^ is the (fc — l)th derivative of the polynomial initial data fa{u)- The function \E'^(r; u) 

in (|4.2|) satisfies a similar linear over-determined system. 



Theorem 4.2 (First Main Theorem) Let he f{u) a smooth monotone increasing function with 
domain {a,b), —00 < a < b < +00. If qk — qkiui,U2, ■ ■ ■ ,M2g+i)j ^ l£ k < g, is the symmetric solution 
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of the linear over- determined system 

^d^qk{u) dqk{u) dqkiu) 



2{ui - Uj) 



duiduj 



du. 



i^j,i,j^l,...,2g+l, g>0 



qk{u,u, ...,u) = Fk{u) 

2g+l 



(4.10) 



Fk{u) 



2(9- 



(25-1)!! du3~^ 

with the ordering ui > U2 > ■ ■ ■ > U2g+i, then 



1 



2g+l g k 



k=l 



1=1 



i = l,...,2g + l, 



(4.11) 



solves the boundary value problem {3^}, {3.1 - 3. It ). Conversely every solution of (3^4), {3/1- 3.1 d^) can 
be obtained in this way. 

Before proving the theorem we show how to obtain the solution of system ( |4.10 ) for generic smooth 
initial data. We follow the procedure in |^ . We start with the following lemma. 

Lemma 4.3 ^ The syste 



tem 



2{z-y)p^y = Pz-pPy, P>0 
p{z,z) = s{z) 



(4.12) 



has, for any smooth initial data s{z) one and only one solution. Moreover, the solution can be written 
explicitly 



P{z,v) = 



where 



P J-i 



(4.13) 



(4.14) 



Using the above lemma, the linear over-determined systems (4.10) can be integrated for any smooth 
initial data in the following way. Suppose that qk{ui, U2, ■ . ■ , U2g+i) is a solution of (4.10). 
Clearly ^fe(-ui,M2g+i) = qk{ui,ui, . . . ,ui,U2g+i) satisfies 



2g 



2{ui - U23+1) 



duidu2g+i dui 
Ak{u,u) = Fk{u) 



dAk ^ dAk 



du 



2g+l 



(4.15) 
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which by lemma |4.3| imphes that 

Ak{ui,U2g+l) 



1 FkC-±I^Ul+'-^U2g+l) 



C^g 7-1 y/1 - £,2g 

For each fixed W2g+i the function Bk{ui, U2g, U2g+i) — qk{ui, ■ • • , wi, U2g, U2g+i) satisfies 



(4.16) 



2{UI - U2g) 



Bk {u,U,U2g+l) 

Using again lemma we obtain 



23-1 

d^Bk ^ dBk dBk 
duidu2g dui du2g 

= Ak(u,U2g+l) 



(4.17) 



Bk{ui,U2g,U2g+l) / / d6grf6s-l(l+6g)^ '(l+6ff-l)' 



^k[^—( 2 "1 ^ 2 U2g)^—U2g+l) 



2 X 



(4.18) 



Going on in the process of integration we obtain the solution qk{u) = qk{ui,U2^ ■ ■ ■ ,U2g+i) of the 



boundary value problem ( 4.10 ) namely 

r-l rl pi 



Qk{u) =i J ^ j ■ ■ ■ J ^ dCid6 . . . + ^2gr-\l + ^2g-iy'^ • ■ • (1 + 6)^ (1 + 6)"^ 



(4.19) 



Fk{^{. . . {^{^U, + i^,i2) + ^U3) + ...)+ '-I^U2g+l) 



V(i~ei)(i-6)-.-(i-6s) 



where C = Y[j=i has been defined in ( 4.14 ). When the initial data is of the form ( 3.11 ), 

the expression (gl|) for the (7fc(u)'s is equivalent to ( [l.2| ). 

Theorem 4.4 The boundary value problem ( 4-l(\ ) has one and only one solution. The solution is 
symmetric and is given by ( ^.Ij ). 



Proof: Uniqueness follows from lemma 4.3 and the argument previous to (4.19). The boundary 
condition (4.10) is clearly satisfied. The symmetry follows from the construction. Indeed in the process 
of integration we can interchange the role of any of the variable Ui . □ 

We have the following relations. 
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Lemma 4.5 The functions Fk{u) and the solutions qkiu), k = I, . . . , g, of the boundary value problem 

g relations. 

-Fk+i{u) + uduFk+i{u), k^l,...,g-l, g>0 



{4-. It ) satisfy the following relations. 

2g + l 



duFk{u) 



(4.20) 



1 



duiQkiu) = ^qk+i{u) + u,duiqk+i{u) i = 1. . . . , 2g + 1, k^l,...,g-l g > 0. 



Proof of Theorem |4.2| (First Main Theorem). 

We consider the non trivial case where qk{u) ^0, k = 1, . . . g, and du^Qgiu) ^ 0, j — 1, . . . ,2g + 1. 

The proof consists of three parts. 

a) The Wi{ii) 's defined in ( j-ll ) satisfy l[3.^ ). 

Using the definition of Wi{u) in (4.11) we have the following relation 



dujWi{u) =2- 



' a 



-dujdu,qg{u) - 2- 



Pi in,) 



-du,qg{u) 



-du, [ ^(2Z ~ 1) ^-/^^^j^ ^gfc(77)rfc_; ) , i^j, *,j = l,..., 2.9 + 1. 



(4.21) 



k=l 



The following identities hold g] 



duj PSiUi) X^ - Xj dUj \Pg{Ui) Pi{uj) 



i^j,ij = l,...,2g+l, k>l 



(4.22) 



1 9A, du^Piiu,) 1 1 



\i — \j duj 



Pi{u^) 



2 Uj — U4 



* 7^ j = 1, • • ■ ,2g+ 1, 



(4.23) 



where Xi{u) has been defined in (3.1) and P^{r) has been defined in (2.4) 



Using the definition of Tk in ( 2.19 ) it is easy to verify that 



dTk 

duj 



k-l 



(4.24) 



Applying repeatedly the relations ( 4.20| ) we obtain the following expression for du^qkiu). 



g-k 



1 \. 

du,qk{u) = ^ X! Qrn+k{u)uf^'^ + u^f''du^qg{u), k = l,...,g-l. 



(4.25) 
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From (4.10) and (4.22- 14.251 ) we can write dujWi{u) in (4.21) in the form 
]fliX]^iiui - ui) du,qg{u) - du^qgiu) 



dujWi (u) 



-2 :5r^ — ^ Ou,qg[u) 



E(2;-i) 

1=1 

E(2^ - 1) 



Pl-IM / 1 



Pq{u%) Ui-Uj 
9-1 g-k 



P(1M \2 



k—l ni—1 



qm+k[U)Uj 
k-l 



k=l 



1=1 
g 



P^M 



k=l 



^ / - k — l — m^j 

2 ^ 

m=l 



1 dX, (PlM) Pf_,{u,) 



1=1 



Simplifying we obtain 



E 9fc(M)ffc-i i 7^ j, «, j = 1, . . . , 25 + 1. 



(m, -Uj)Po^(ui) 



dujqgiu) + 2- 



dmqgii 



1 aA, 

Ai — A, Sm, 



1=1 



1 

Ai — A, i9u 



E(2'-i) 



Pi-iM 



E qk{u)fk-i j, i,j = 1,..., 2.9 + 1. 

(4.26) 



k= 



Adding and subtracting the quantity — — — Wj to (4.26), we can reduce it to the form 

Ai — Xj ouj ' ^ 

du,W^-- - — [W^ - W,] = du,qg{u) { - — ^ 



' Xi — Xj duj 

T23+1 



Xi — Xj duj P^{uj) 



Y{l^l^^^{u^ ~ Ul) 
{Ut-Uj)P^{ui) 



(4.27) 



k-iw 



g-k 



k=l 



The term in parenthesis in the right hand side of ( 4.27 ) does not depend on the initial data f{i 



It is identically zero for the analytic initial data (3.11) because in such case the WiS satisfy (3.4) 
Therefore we can conclude that 

1 dX 



Ai - A, oui 



(4.28) 



for any smooth monotone increasing initial data x = f{u). 
h)The Wi{u) 's satisfy the boundary conditions (3.'/- \3.ld^ . 
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In the following we use the sup-script g to denote the corresponding genus of the quantities we are 
referring to. We have the following relations. When ui = u;+i = v for 1^1,..., 2g, the F^'s defined 
in ( |2.19| ) satisfy 



f^(ui, . . . ,Ul-i,V,V,Ul+2, . . ■,U2g+l) =ff ^(Ui, . . . ,U/_i,M;+2, ■ • • ,U2g+l) 



vtl_\{ui, . . . ,U/-l,U;+2, . . . ,U2g+l), fc > 1, 5 > 1 



(4.29) 



and the (7fe(-it)'s defined in (4.19) satisfy 



qliui, . . . ,Ul-i,V,V,Ul+2, ■ ■ ■,U2g+l) - Vql^^{ui, . . . ,Ul-i,V,V,Ui+2, ■ ■ • ,U2g+l) 



= ql ^(ui, . . . ,M;_i,w/+2, ■ • ■ ,U2g+i), k ^ 1, . . . , g - 1, g > I. 

(4.30) 



For Ui =/= ui ^ ui^i = V we have that 



dinQg-liui, ■ ■ . ,U/_i,U;+2, • ■ ■ ,U2g+l)--q'^{ui, . . . , Ui_i, V, V, Ui+2, ■ ■ • ,W2g+l) = 

{Ui - V)duiq^g{ui, . . . ,Ui^i,V,V,Ui+2, ■ ■ ■,U2g+l), 



(4.31) 



which follows from (4.20). When ui = = v the polynomials P^{r)^s defined in (^J) satisfy the 
relation |19l 



P^{r) = {r -v)P^^^{r), k>0. 



(4.32) 



Using the relations (4.2E-4.32) we have that for z 7^ ^, Z + 1, z = 1, . . . , 2g + 1, 

n2g+i / _ \ 



wf(ui, . . . ,Ui_i,V,V,Ul+2, ■ ■ ■,U2g+l) = 2 



m=l U V k=m 

/ 9 



-du^Tg-l 



5-1 



Vm— 1 



which reduces to the form 



(Ui, . . . , Ui_i,V, V, Ui+2, ■ • • , "29+1) 
/ 9 



„9-l 



(Ul, . . . ,U/_i,M;+2, . . ■,U2g+l) 



\^m— 1 



(4.33) 



(4.34) 



Using ( 2.21 ) the term in parenthesis in the right hand side of ( 1.34 ) turns out to be identically zero. 
Therefore the boundary conditions (p^) are satisfied for any smooth monotone increasing initial data. 
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Since the functions qk{uys in ( 4.ig| ) and the ffe({t)'s defined in ( 2. 19] ) are symmetric with respect to 
Ml, M2, . . . , U2g+i we immediately deduce from ( [4.32[ ) that 



wf (Ui,. . . ,Ui^i,V,V,Ul+2, ■ ■ -,"25+1) = wf+i(ui, . . . ,M;_l,W,U,U/+2, • • ■ ,U2g+l) 



(4.35) 



so that the boundary conditions (3.7) are satisfied (for a more detailed analysis of this limit see 
section 5^). When 5 = 1 we deduce from ( 4.35 ) 

Wl(ui,Ul,U3) = W2(m1,Ui,U3) 



and from ([4.3lH4.32|) 

W3(ui,Ml,U3) = 2(U3 - Ui)l9„3gi(Mi,Mi,M3) + gi (ui, Ui, M3). 

From (|4.10D and (|4.19D we get the relation 

gi(ui,ui,U3) = /(U3) + 2{ui - U3)5„3 9i(ui,ui,U3) 

so that 

W3{ui,Ui,U3) = fius). 



An analogous result can be obtain when U2 = U3, so that the boundary conditions (3.9-3.10) are 
satisfied. 

c) Uniqueness. We prove that when f{u) = then wf (u) = for all > ui > 1*2 > • ■ ■ > U2g+i > a, 
for 1 < i < 2g + 1 and for any 5 > 0. The proof is by induction on g. For g — the statement is 
satisfied. 



For 5 = 1 we repeat the arguments of ||6|. We fix U2 and we consider the equation (3^) with 
boundary condition (3.9|-3.1C), namely 



du3 



ai3[wi - W3] 



dw3 r 1 

= a3i[W3-Wi] 

OUi 

Wi(ui,U2,U2) = /(ui) = 

W3{U2,U2,U3) = /(m3)eeO. 

We can regard the above equations as a first-order linear ordinary differential equations with non 
homogeneous term. Integrating them we obtain a couple integral equation. By standard contraction 
mapping method it can be shown that when /(u) = this system has only zero solution, i.e. wi — 
W3 = for (iti,M3) satisfying > ui > U2 > ^3 > a. Because of the arbitrariness of U2, wi and ^3 
vanish as a function of (ui, U2, U3) and therefore, by (|3.4[) so does W2{u). Now we suppose the theorem 
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true for genus g — 1 and we proof it for genus g. We fix 6 > U2 > U3 > • ■ ■ > U2g > a and we consider 



the equation ( |3.4[ ) for wf and W2g^i with boundary condition (3.7-3.8), namely 

d 



n9 — 



dU2g + l 

d „ 



ai(2g+l)K - iwfg+l 



J(2g+l)l[W2g+l 



wf{ui,U2, ■ . .,U2g,U2g) 
W2g+l{u2,U2,...,U2g,U2g+l) = U'fg;^i(u2,U2,U3, ■ • • ,W2g,M2g+l) 



(4.36) 



; 
0. 

= w: 



2g+ 



^{u) = 0, for 



Repeating the arguments developed for genus 17 = 1 we may conclude that wf (u) 
arbitrary b > ui > U2 > ■ ■ ■ > U2g+i > a- We then repeat the above argument fixing 6 > ui > U3 > 
• • • > U2g-i > U2g+i > a and considering the equations (3.4) for (w) and W2g{u), namely 

d 



dU2g 

du2 



w?, 



"2s 



«2(2g) 



W?, 



W: 



2g\ 



W^{ui,U2, . . . ,U2g-l,U2g+l,U2g+l) 
W2giui,Ul,U3, . . . ,U2g,U2g + l) 

It can be easily shown that also ^^(u) 



a(2g)2[«-'2g " W^] 



(4.37) 



Ujf (wi,U2, . . . ,U2g^l,U2g+l,U2g+l) = 



"fg ^{U1,U1,U3, . . .,U2g,U2g+l) = 0. 

for arbitrary 6 > ui > U2 > • • • > "29+1 > a- 



W2g{u) 

Repeating these arguments other g — 2 times we conclude that wf (u) = ior 1 < i < g, g + 2 < i < 



2g + 1 and for arbitrary b > ui > U2 > ■ ■ ■ > U2g+i > a- Applying (B.4) and the boundary conditions 
(3.7)-(3.8) we can prove that also Wg^i{u) is identically zero. The theorem is then proved. □ 

In the next section we consider the problem of reality of the solution of the hodograph transform 



(3.6). 



Remark 4.6 The boundary conditions (S.1- 3.1 guarantee the -smoothness of the solution of the 
Whitham equations. Indeed it can be proved that the x derivatives dxUi{x,t), i = I, . . .2g + I, are 
continuous on the phase transition boundaries. 



5 An upper bound to the genus of the solution 

In this section we give an upper bound to the genus of the solution of the Whitham equations for the 
monotone increasing smooth initial data x = /(u)|t=o- 

Theorem 5.1 (Second Main Theorem) If the monotone increasing smooth initial data f{u) sat- 
isfies the condition 



jf{u) > 0, 1 < iv e IN, 



(5.1) 



for all u e (a, b) except at one point, then the solution of the Whitham equations ( LJ_ ) has genus at 
most N . 



19 



Proof: the solution of the Whitham equations (1.1 ) for different .g > determines a decomposition of 
the x~t plane, i > 0, into a number of domains Dg with g=0,l,2 . . . , (see Fig. 5.3). To the inner part 




X 



Figure 5.3: An example of decomposition of the x — t plane 
of each domain Dg it corresponds the g-phase solution b > ui{x,t) > 1*2(2;, t) > • • • > U2g+i{x,t) > a 



of the Whitham equations (1.1). The common boundaries of the domains Dg, g > 0, and I?g+„, n > 1, 
are the points of phase transition between the g-phase solution and the {g + n)-phase solution. 

We will show that each domain Dg, g < N, does not have a common boundary with any of the 
domains Dm, m > N. Since the set of domains {Dg}g<N is not empty because ^ the set 
{Dg}g^]\f must be empty. Indeed on the contrary the x ~ t plane, t > 0, which is a connected set, 
would be split into a number of domains whose union forms a disconnected set. 

Before determining the boundaries of the domains Dg, g > 0, we need to study more in detail the 



hodograph transform (3.6). 

Proposition 5.2 Let us consider the polynomial 

ZS{r) := -xP^{r) ~ 12tPf{r) + RS{r), (5.2) 
where _R°(r) = /(u) and R^{r), g > 0, is given by the expression 

2g+l 2g+l g k 

R^{r)^2Y,d^,qg{u) n {r-Ui) + Y,qt,{u)Y,{2l-l)f^_iPl,{r), (5.3) 

k=l l=lj=ik k=l 1 = 1 



with the polynomials Pj'{r), I > 0, defined in (JSj) and the functions qk{u), k — 1, . . . ,g, defined ir 



( j.lL ). Then the hodograph transform (3A_ ) is equivalent, for g > 0, to the equation 



Z3{r) =0, g> 0. (5.4) 
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Proof: We observe that the Wi{u)^s defined in (4.11) are given by the ratio Wi{u) 



l,...,2g + 1, where R^{r) is the polynomial defined in (5.£). Hence we can write the hodograph 
transform ( ) in the form 

[-xP^{r)-l2tP({r)+Ra{r)]r^,,^^Q, i ^ I, . . . ,2g + I. (5.5) 

For g > 0, Z3{r) is a polynomial of degree 2g and because of (|5.5D it must have at least 2g + 1 real 



zeros. Therefore it is identically zero. Hence for g > 0, ( p.4[ ) is equivalent to (5_^) and ( |3.q ). □ 

In the following analysis we give some conditions for the existence of a real solution b > ui{x, t) > 
U2{x,t) > ... ,U2g+i{x,t) > a of the hodograph transform (3.6). For the purpose let us consider the 
function '^^(r; u) which solves the boundary value problem 



d 



\E'9(r; u) 



d 92 

*9(r;w) = 2(u, -Uj)^— ^*9(r;u), ij^j, ij = l,...2g + l 



du 



du 

"a^ir^u) - 2^— ^'f(r;u) = 2(r - -u-,) ^^^ ^g(r;M), j = 1, . . .2^ + 1 



drdui 



(5.6) 



^'f(r;r, . . . ,r) = 



29 



2g+l 



(25+1)!! 



where f''^\r) is the f;th derivative of the smooth monotone increasing initial data f{u). From the 
results of Sec. |4| we are able to integrate (5.6) obtaining 



"f'ir; u) =1 J'^ j^- ■ ■ J'^ . . . rf63+i(i + C2<,+i)^(i + ^23)'-^ ... (1 + 6 



2 X 



(5.7) 



2 



U2) + ...) + 



— W2s+l) 



v/(i-a)(i-6)...(i-69+i) 



where X 



29 



7 rrr TTw^t^ C*?" and the C^'s have been defined in (kt.l4|). The function '^^(r:u) is 

{2g + l)\\'-'-^=^ ^ ' ^[^^ ^ ' ' 

symmetric with respect to the variables 5 > mi > U2 > ■ ■ ■ > fJ-2g+i > a. For the initial data (3.11), 
the expression of ^'9(r;u) in (5.7) is equivalent to (4.2). 

Proposition 5.3 Let us consider the function 



2g+l 



(5.8) 



where 'ii^{r;u) has been defined in (5.1). Ifui{x,t) > U2{x,t) > ■ ■■ > U2g+i{x,t) satisfy the Whitham 
equations (1.1) then the function ^^{r;u) has, in the r variable, at least one real zero in each of the 
intervals {u2k, U2k-i), k — 1, . . . , g, g > 0. 
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Proof: If the functions ui(x,t) > U2{x,t) > • • • > U2g+i{x,t) satisfy the Whitham equations then 
by proposition 5.2 the polynomial Z^{r) = 0. Therefore 



Z9(r) 



dr = 



-xP^{r) - 12tPf (r) + Ra{r) 



dr 



/i(r) 



(5.9) 



In the third equality of (5.£) we have used the fact that 



dr 



crf(r) =0, l>0, k = l, 



,9, 



because of the normalization conditions (2.5). The function '^^{r; u) satisfies the relations 



r - Ui 



(5.10) 



which can be easily obtained from ( ^.7| ). Using ( ^.1C| ) we can rewrite the last term in (|]^) in the form 



=2 



/i(r) 

(r; u) 
r - Ui 



-dr. 



2a„,*s(r;w) 



^l{r) 5,*s(r;?I) + £9„.^'f(r;u) dr, fc = l,...,g, 



i=l / 

where the last equality has been obtained integrating by parts. Using the definition of (f>3(r; u) in ( |3.8[ ) 
we rewrite the above relation in the form 







fi{r)^^{r;u)dr, k = l,...,g. 



(5.11) 



Relation (5.11) is satisfied only if the function <I>3(r;u) changes sign at least once in each of the 
intervals (u2fe,W2fc-i), k ^ 1,.. . ,g. □ 

In the following we are going to determine the equations which describe, on the a; — t > plane, 
the boundary between the domains Dg and 

The boundary between the domains Dg and -Dg+i represents a singular behavior in the solution of 
the g or (g+ l)-phase equations. As we have shown in the example on fig LI , the boundary between the 
domains Dq and Di is described by the curves x^{t) where U2{x,t) = U3{x,t) and ui{x,t) = U2(x,t) 
respectively and by the point xq , to of gradient catastrophe of the zero phase solution. 
In a similar way the generic boundary between the domains Dg and -Dg+i is described by 
a) the curves x^{t) where two Riemann invariants of the (g + l)-phase solution coalesce; 
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b) the point of gradient catastrophe Xc, > of the g-phase sohition, namely the point where one 
of the 2(7+1 Riemann invariants has a vertical inflection point. The equations determining the point 
of gradient catastrophe of the g-phase solution can also be obtained considering the limit of the g + I 
phase solution when three Riemann invariants coalesce. 

To treat case a) we consider the Riemann surface iSg+i of genus g + 1 given by the equations 

jl'^ ^ {r - V - VS){r - V + y/S)ii^, wGlR, 



2s+l 

A*^ = - itj), b > ui > U2 > ■ ■ ■ > U2g+i > a, 

i=i 

where v ^ Uj, j = 1, . . . 2g + 1, and < (5 <C 1. The Riemann invariants are the 2g + 3 variables 
ill — V + \/6, U2 = V — \/6, ui > U2 > ■ ■ ■ > U2g+i- We suppose ui,U2 ^ Uj, j — 1, . . . ,2g + 1. 
The hodograph transform ( |3.6| ) for these 2g + 3 variables has two different behavior for S when 
V belongs to one of the bands 



{u2g + l,U2g) U iu2g~l,U2g^2) U • • • U (u3,M2) U 



(5.12) 



or gaps 



(a, U2g+l) U {U2g, U2g-l) U • • • U (W4, Us) U (U2, Ui). 



(5.13) 



We call leading edge of the phase transition boundary the case in which v belongs to the bands. We 
call trailing edge of the phase transition boundary the case in which v belongs to the gaps ( 5.13| ). 

Theorem 5.4 The leading edge of the phase transition boundary between the g-phase solution and 
the {g + l)-phase solution is described by the system 



$9(„;t/) _ Qtego = 
a^$f(w;M) = 

P^{r) ^ P^{r) 



-I2t 



(5.14) 



1 = 1, 



,25+1, 5>0 



where v € (w2j+i,U2j), < j < g, uq = +b, the function ego = 1 for g = and zero otherwise, the 
function $^(r;'u) has been defined in ^5.^ and the polynomial R3{r) has been defined in ^5.^ . We 
assume (dy)'^^^ {v; H) ^ and ^^{ui; il) ^ 0, i — 1, . . . ,2g + 1, on the solution of (5.14 )■ 



Remark 5.5 System (5.14) is a system of2g + 3 equations in2g + A unknowns x,t,v and ui> U2> 
■ ■ ■ > U2g+i. If system (5.14) is uniquely solvable for real x, v and ui > U2 > ■ ■ ■ > U2g+i as a function 
oft>0, then a phase transition between the g-phase solution and the [g + \)-phase solution occurs. 
The curve x~ — x'Z {t) describes on the x~t > plane the boundary between the domains Dg and -Dg+i 
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associated to the leading edge. The conditions u) ^ and ^^{ui\ u) 7^ 0, i = 1, . . . , 2g + 1, 

on the solution of (5.14) exclude higher order degeneracy in the transition. Indeed in such case it can 
be proved that in a neighborhood of the solution x{t), v{t), ui{t) > U2{t) > • • • > U2g+i{t) of ( 5. 1^ ) 
the {g + l)-phase solution is uniquely defined. 



Theorem 5.6 The trailing edge of the phase transition boundary between the g -phase solution and 
the {g + l)-phase solution is described by the solution of the system 

$f(t;;u) - 6t€go = 

($3(r;ii) - 6tego)n{r)dr = 
Pf(r) R3{r) 



(5.15) 



-t 



P^ir) P^ir) 



i = l,...,2g + l, g>0 



where v £ {u2j, U2j-i), 1 < J < 5 + 1, U2g+2 = a and the function $3(r; u) has been defined in (5.i ). 
We assume w) 7^ and ^^{m; m) 7^ 0, i = 1, . . . , 2^ + 1, on the solution of ( 5.15 ). 

If system ( ^.15 ) is uniquely solvable for real x, v and ui > U2 > ■ ■ ■ > U2g+i a-s a function of t > 0, 
then a phase transition between the g-phase solution and the (g+ l)-phase solution occurs. The curve 
x^ = x^{t) describes on the x — t > plane the boundary between the domains Dg and -Dg+i 
associated to the trailing edge. 

The following theorem enables one to determine a point of gradient catastrophe of the g-phase solution. 
This point can be obtained either as a limit of the (.g+1) phase solution when three Riemann invariants 
coalesce or imposing a vertical inflection point on the g-phase solution. 

Theorem 5.7 Let us consider the system 



dr<i>a{r;u)\r=u, =0 
<i>s{ui;u)-6tego^O 



1,...,25 + 1, 5>0 



P^ir) Pi{r) 

where CgQ = 1 for g — and zero otherwise, the function ^^{r;u) has been defined in ^ 



(5.16) 



When 



system (5. It) is uniquely solvable for Xc, tc > and b > Ui{xc, tc) > U2{xc, tc) > ■ ■ ■ > U2g+i{xc, tc) > 
a, then the g-phase solution has a point of gradient catastrophe on the ui branch, 1 < I < 2g + 1. 



To avoid higher order degeneracies in the transition we impose the condition 

(a,)2$9(r;w(a;„t,))|,=„,(,^,i^) 7^0. 



(5.17) 
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Indeed we can prove that the above condition guarantees that the genus of solution of the Whitham 
equations increases at most by one in the neighborhood of the point of gradient catastrophe. Therefore 
it is legitimate to consider the point of gradient catastrophe that solve ( 5.16| ) and satisfies ( 5.17 ) as 
a point of the boundary between the domains Dg and i'g+i. The condition (5.17) in not essential in 



the genus g — case as illustrated in Theorem 3.2 



Remark 5.8 We observe that both systems (5.14) ( 5.l9( ) in the limit v —>■ ui, I = 1, . . . ,2g + 1, 
coincide with system (5. It ). 



Theorems 5.4, 5.6 and 5.7 characterize all types of boundaries between the g-phase solution and the 
(3+l)-phase solution, namely the leading edge, the trailing edge and the points of gradient catastrophe. 
We will prove these theorems in the next section. 

Example 3.1 For x — k — 3, 5, 7, ... , the solution of the Whitham equations has genus at 
most equal to one |l^,|6j. On the x — t, t > 0, plane we have only the domains Dq and Di. The one 
phase solution is defined within the cusp a;_(fc) < x < Xj^{k) tT^ , where x_(k) < Xjf.{k) are two 
real constants and t > 0. The point x = 0, i = is the point of gradient catastrophe of the zero-phase 
solution. The constants a;_(fc) and x^{k) can be computed explicitly. From (5.7) and (5.8) we obtain 



$o(?', u) 



1 



On the leading edge where ui = U2 = v and — v > ( ^.14 ) is equivalent to the equations 



(5.18) 



/ 

J U 



no - 6t 



d^^O, x = -6tu + f{u). 



(5.19) 



On the trailing edge where ui = u, U2 = = v, u > v, system ( ^.15 ) is equivalent, for g = to 

{f'{0-(it)V^d^^Q, I iM^d^^Q, x = -6tu + f{u). (5.20) 

Equations ( 5.19 ) and ( ^.2C| ) have already been obtained in [||. Solving system ( 5.19 ) for t when 
f{u) = ^t^ fc = 3,5,7, ... we obtain x_(/c) (l3|,0 



x_{k) = -6^— ^(2z_(fc) - 1) 
k 



6 



(l + 2(fc- l)2_(fc)) 



(5.21) 



where z-{k) > 1 is the unique real solution of F{—k + 2, 2, |; z) — 0. Here F{a, b, c; z) is the hyperge- 
ometric series. The quantity X-^-{k) is obtained from (5.20) 



fc — 1 

x+(k) = 2——(2z+(k) - 3) 
k 



-(3 + 2(fc-l)z+(fc)) 



, fc — 3, 5, 7, . 



(5.22) 
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where the number z+(fc) > 1 is the unique real solution of the equation F{—k + 2. 2, |; z) = 0. We 
give some numerical values: a;_(3) = — 12\/3, x-{5) = —16.85, x^{7) = —16.21, x-{9) = —16.09, 



x+(3) 



a;+(5) = 1.58, x+{7) = 1.61, x+i9) ^ 1.72. 



Theorems 5.4, 5.6 and 5.7 can be generalized to the case of multiple phase transitions. 

Theorem 5.9 The transition boundary between the domains Dg, g > and Dg+n, n > 1, having 
Hi leading edges and n2 trailing edges and points of gradient catastrophes, ni + n2 + ns = n, 
ni,n2,n3 > 0, is described by the solution of the system 

c»t,^^$9(i;fc;u) = 0, k^l,...,ni 
<i>^{vk;u)-6tego = 0, fc=l,...,ni 
<^>s{yi;u)-6t€go^0, l = l,...,n2 
C"~"i'^'ir;u) - 6tego)tiir)dr = 0, 1^1,. 
dr^3{r;u)\r=uj^ = 0, to = 1,..., 713 
<i>^ (wjm ; u) — 6i Cgo = m = 1 , . . . , na 



,n2 



(5.23) 



Pf{r) R3{r) 



2= 1,..., 25+1, 



where Vk & (^2^^ , M2i;,+i); < jk < g, k ^ 1, . . . ,ni and yi e {u2j,-i,U2j,), I < ji < g + I, 

Z = 1, ... ,712. 



Remark 5.10 The transition between the g-phase solution and the {g + n)— phase solution, n > 1, 
is highly non generic (cfr. Fig 5.2o ). Indeed system (5.2.1 ) is a systems of 2n + 2g + 1 equations in 
ni + n2 + 2g + 3 unknowns Vi,.. . ,Wni, • • ■ ,2/ri2 7 "i, . . . ,1*23+1, x and t > 0. Therefore if system 
(5.2 J ) admits a real solution, n + — 2 variables are functions of all the others. 



Theorem 5.£ includes all the degenerate cases that have been excluded in Theorems 5.4, 5.6 and 5.7. 
For example let us consider a transition with a double-leading edge described by the variables vi, V2 
and ui > U2 > ■ ■ ■ > U2g+i which solve ( 5. 2^ with ui = 2, n2 = 0, = 0. In the limit vi ^ V2 — v 
we get a .single degenerate leading edge. The variables v , ui > U2 > ■ ■ ■ > "29+1 satisfy ( |5.4 J o,nd the 
equations (dj,)^^^ {v;u) = and {dv)^^^{v;u) = 0. Therefore we regard such degenerate single leading 
edge as a point of the boundary of the domains Dg and Dg+2 ■ 



Proposition 5.11 // the domains Dg and -Dg+i have a common boundary, then the function 

^3{r;u)-&tego 

has at least g + 2 real zeros in the r variable for Ui{x, t) > U2{x,t) > ■ ■ ■ > U2g+i{x, t) satisfying l[5.^ 
and for t > 0. 
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Figure 5.4: an example of phase transition between the one-phase solution and the 4-phase solution 
having one leading edge, one trailing edge and one point of gradient catastrophe at U3. 



Proof: The domains Dg and the -Dg+i have a common boundary if one of the systems ( ^3.14 ), ( 5.15 ) 
or (5.16) is solvable for some t > 0. We first consider the leading edge. For g — Q the statement 



is obvious from system (5.14). For g > system (5.14) can have a solution if the function $^(r;u) 
has a double zero at r = v when v belongs to the bands ( 5.12 ). Combining this observation with 



proposition 5.3 we immediately obtain the statement. As regarding the trailing edge, the theorem is 



obvious for g = 0. For g > and v e (a, U2g+i) system ( 5.15| ) can be satisfied if the function <i>s(r; u) 



has at least two zeros in the interval (a,U2g+i). Combining this observation with proposition 5.3 



we 



obtain the statement. When v € {u2j,U2j-i), 1 < j < 5, (5.11) and (5.15) are satisfied if the function 
^s{r;u) has at least three real zeros in the interval (u2j, W2j-i), 1 < i < ff, and changes sign at least 
once in each of the intervals {u2k, U2k-i), k — 1, . . . , g, and k ^ j. Therefore ^^{r; u) has at least g + 2 
real zeros at the trailing edge. If the point of the boundary between the domains Dg and the -Dg+i 
corresponds to a point of gradient catastrophe of the g-phase solution the statement is obvious from 



system ( 5.16 ) and proposition 5^. Therefore, on the phase transition boundary between the domains 



Dg and the function ^s(r;u) has at least g + 2 real zeros in the r variable, b > r > a, 

b > Uiix, t) > U2{x,t) > ■ ■ ■ > U2g+i{x, t) > a and for < > 0. . 



for 

□ 



Proposition 5.11 can be generalized to multiple phase transitions. 



Proposition 5.12 On the phase transition boundary between the domains Dg and the Dg-^-n the 
function 

$f(r;u)-6tego 

has at least g + 2n real zeros in the r variable for b > ui{x,t) > U2{x,t) > • • • > U2g+i{x,t) > a and 
for t > 0. 



The proof is analogous to that of proposition 5.11 
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Lemma 5.13 If the smooth initial data x = /(u)|t=o satisfies ( |5. then the function 

has at most 2N — g real zeros (counting multiplicity) in the r variable for < g < 2N and for any 
real x, t > and b > ui > U2 > ■ ■ ■ > U2g+i > a. 

Proof: for proving the lemma we need the following elementary result. If the real smooth function 
^(r) satisfies the condition 

— ^(r)>0, 0<me]N, 
for all r belonging to the domain of ^, then ^(r) has at most m real zero. 



Using ( p.l\ ), ( |5.6| ) and (^) we obtain 

^^'^'ir- u) = j ■■■ j ^di,... (1 + C2,+i)^^(i + 6.)^^-^ • . . (1 + 6)^^-''+^ (1 + ^if-' 

2(2s+l)(2iV-s) K V(l-6)(l-e2)...(l-6s+l) 

(5.24) 



for (? > 0, for any real r and for any fixed real iti > U2 > • • • > W2g+i belonging to the interval (a, b). 
Therefore $3(r; u) has at most 2N — g real zeros in the r variable for < g < 2N and for any fixed 
real 6 > wi > U2 > ■ • • > M2n+i > a- D 



Lemma 5.14 If the smooth monotone increasing initial data x — f{u)\t=Q satisfies (5.1), then the 
function 

$"(r;M) - 6t 

has at most 2N real zeros (counting multiplicity) in the r variable, for any t > and u ^ u* , where 
u* is the unique solution of the equation f^^^^^'>(u*) = 0. For what ^^{r;u*) — 6t is concerned we 
have two possibilities: either it has at most 2N real zeros (counting multiplicity) in the r variable for 
any t > and some of these zeros are distinct, or it has at most two real zeros in the r variable for 
any t > except for t = t* > where r = u* is a zero of multiplicity higher than 2N . 



Proof: Since f'{u) > and because of (5.1) we have two possibilities: either f'{u) has at most N 
simple minima and some of these minima are distinct, or f'{u) is a non negative convex function with 
a minimum at m = u*, where u* is the unique solution of the equation /(^^+^)(m*) = 0. From the 
above considerations the lemma easily follows. □ 
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We remark that when f'{u) is a convex function, we can apply Tian's Theorem |3.2| . Indeed in such 
case we have f"'{u) > for aU u ^ u* and therefore the solution of the Whitham equations has genus 
at most one. 

We continue the proof of the second main theorem. In the following we exclude the case in which 



f'{u) is a convex function. From proposition 5.12 we deduce that when the domains Dg and D^, 
> 9 ^ 0, have a common boundary the function <I>^ (r ; u) — Gtego defined in ( ^.^ ) has at least 2m — g 



real zeros for some values of x and t. From lemma 5.13 the function $^(r;u) — 6tego defined in ( 



has at most 2N ~ g real zeros for g < 2N when the initial data f{u) satisfies (5.1). Therefore 



2N ~ g>2m~ g or m< N. 



(5.25) 



This shows that the set of domains {13g}o<g<A' does not have common boundaries with any of the 
domains in the set {D„i}m>N- The second main theorem is proved. □ 



5.1 Phase transitions 



In this subsection we prove theorems 5^, 5^, 5/7 and (5.£). 

In the following wc denote with a sup-script the genus of the quantity we are referring to. Whenever 
we omit the sup-script we are referring to genus g quantity. 

We denote with cr^^^ = a^~^^ {r, S^v) the normalized Abelian differential of the second kind with 
pole at infinity of order 2fc + 2 defined on the surface Sg+i of genus g + 1 



fl^ — (r — V — VS){r — V + v^)/i'^, 



23+1 

= (r - Uj), Ui> U2> ■■■> M2ff+1, 

where v ^ Uj, j = 1, . . . 2g -|- 1, and < 5 <C 1. Wc define the polynomials 

Pk (''' = ■ 

The expression of the polynomial in ( |5.3| ) for the Ricmann invariants ui = v + U2 ^ V — VS and 
ui > U2 > • • • > M2g+i reads 

2g+l 2g+l 

RS+'ir, S,v)^2j2 duAlXl n - ":')(^ - - "2) + 2^l\r){r - u^)du,qlX\ 

+2M^(r)(r - U2)du,qlX\ + ^ q^' Y.^21 - l)f ^,^JPf_Y i^, u, S) 

k=l 1=1 
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where ql^^ = -U2, wi, . . . , U2g+i), k 

We will SOI 
on the parameters v and 5. 



tl^^{ui,U2,Ui, . . . ,U2g+l), 

k > 0. We will sometimes omit the explicit dependence of <7^^^{r, S, v), Pf!^^(r, 6, v) and i?^+^(r, v, S) 



,g+ 1 and ff+i 



5.4 



Proof of Theorem 

We write the {g + l)-phase solution (3^) for the variables iti — v + VS, U2 — v 
U2g+i in the form 







-12t 



ui-u2\ Pt\u^,5) Pt\u,,S) 

Pf+\r, S) Rs+^r, 6) 
Pll+\r,5)^ Pt\r,S) 

Pf+\r,S) R9+\r,S) 



Ul > U2 > ■ ■ ■ > 

Pf+\u2,S) , Ra+Hu2,S)' 



pt\u2,s) pr\u2,s) 



+ 



9+1 /T, 



(5.27) 



Pt\r,5) Pnr,S) 



1,...,25 + 1 



The reason to write the hodograph transform (3.6) in the form (5.27) is that system ( ^.27 ) is non 
degenerate at the phase transition ui ~ U2- 

We show that system (5.27) reduces to system (5.14) when ui — U2 or S = 0. 

It is clear from Theorem 4.2 and proposition 5.2 that the last 2g + 1 equations in ( 5.27 ) in the limit 
-Ul — U2 reduce to those of system (5.14). 

For computing the limit iii — > 112 of the first two equations in (5.27) we first need to study the behavior 
of the differentials a^'^^{r, 5,v), k > 0, when S ^ 0. When v belongs to one of the bands ( 5. IS ), the 
expansion of cr^^^ (r, S, v) for 5 reads 



(5.28) 



where cr^(r) is the normalized differential of the second kind defined on Sg with pole at infinity of order 
2/C + 2, al{v) ^ ^^\r=v and = —. The differential Lu3{r) is the normalized Abelian differential of 
the third kind with poles at the points Q^{v) ~ {v, ±/i(w)) with residue ±1 respectively. The explicit 
expression of uj!^{r) has been given in ( p.l4| ) The differential 0((5^)/(5^ has a pole at r = w of order 4 
and zero residue. 



From (5.28) we can get the expansion of the polynomial Pf!^^{r,S,v) — jj.(r) 



dr ' 



namely 



Pt\r, S, v)^{r~ v)P^{r) + -al{v) f,'{v) - (r - «) £ r^-'N',{v) + 0{S'), 



(5.29) 



fc=i 



where P^{r) has been defined in (2/7) and the A^fc(w)'s have been defined in ( 2.18 ). From ( 5.29| ) we 
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can evaluate the following 



r—v:Ly/S 



0(6). 



(5.30) 



To evaluate the first two equations in (5.27) at the point of phase transition we need also the 
following relations 



duidu 



—qgXl{ui,U2,u) 



^s{v;u), 1 = 1,2, 



(5.31) 



where the function ^>^{r;u) has been defined in (5.7). From ( 4.29| -4.3C), (5.30) and (5.31) we obtain 

B9T\ ^ ^-{v; u) 



^^'«(-^...(.;.)£\2._l)f^,,_.^ 



P'M 



PIM 



(5.32) 



Y^ql[u)Y^i2l^l)VU^^, i = l,2. 



k=l 



1=1 



Piiv) 



From ( 2.21 ), ( 5.10 ) and using the definition of <f>^(r;-it) in ( p.q ) the above relation can be written in 
the form 



RB+\ui,5) 



p!W) ' 



i = 1,2, 



(5.33) 



where R^{r,u) has been defined in (5.3). We need to consider also the quantity 
1 ( R3+\u-^,5) R9+^{u2,5) 



r- ^ (w + ^) „ „ g+1 . 

V^ „g + 1, , /7 ., ^nign2gg+l (^il, "2, W) 

Pq (u + V<5, b) . 



+2 



Ul-U2\P^+\ui,5) Pt\u2,6)^ 

[ [priv+vs,s) - priv^vs,s) ) '^-^ ' 

g,,(.)E(2^-i)r,_,(.) i^— 



5=0 



5=0 



9+1 



1 / P;lY + V^, 5) Pf_Y {v-VS,6)\ 



(5=0 



where we have used ( [4.29H4.30| ) to obtain the right hand side. Using ( ^.2l[ ), ( |5.10| ), ( |5.30D and ( ^.3l[ ) 

the above reduces to the form 



1 R9+^Ui,d) RS+\il2,S) 

ui~U2 \p^+\u^,5) " Pt\n2,5) 



(5.34) 
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From ( 5.3C| ), (5.33) and (5.34), the system (5.27) on the point of phase transition ui — U2 = v 
reads 

d f-mPf{v) + R9{v,u) + 2fj.^{v)<i>9{v;uy 







dv \ Pg(i;,M) 



-I2t 



(5.35) 



P^ir,u) Pl{r,u) 



, i = l,...,2g + l 



From proposition 3.3 the last 2g + 1 equations in (5.35) are equivalent to the condition 

-a;P(f(r) - 12iPf (r) + i?f (r) = , g > 0. 



Therefore system ( 5.35| ) is equivalent to system ( 5.14 ) for 5 > 0. For 5 = substituting (3^) in ( 5.35 ) 
we obtain 



^ dv{-6tu + J{u) + 2{v ~ w)($o(w, u) ~ 6t)) 



< X = -mu + f{u) + 2{v -u){^o{v,u) ~Qt) 



(5.36) 



X = —6tu + f{u) 



where ui = u and f{u) is the initial data. It is clear that (p. 36) is equivalent to (5.14) for g — 



□ 



5.6 



Proof of Theorem 

For getting the equations determining the trailing edge we have to repeat all the above calculation 
in a slightly different way. We write the (g+l)-phase solution for the variables iii = v+\^, U2 = v — \f5, 
ui > U2 > ■ ■ ■ > U2g+i in the form 



Pf^\u2,S) Pf+\ui,6)\ Rs+\u^,6) , Rs+\iL2,d) 







l0g(iii^)2(ui-7i2) 



12t 



Pt\^2,5) P^^\u,,S)l P^^'{ui,5) ■ P^^'{ii2,S) 



12t 



-12t 



Pf+\ r, 5) R3+^r,S) 



ir,S) 



(5 



Pf^\r,5) R3+^{r,S) 



Pt\r,S) ■ Pr\r,S) 



,25 + 1 



The reason to write the hodograph transform in the form ( 5.37 ) is that system ( 5.37 ) is non degenerate 
at the phase transition ui — U2 — v when v G (w2j, U2j-i), 1 < i < .9 + 1, = «• 
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From Theorem 4.2 and proposition 5.2 we deduce that the last 2g+l equations in (5.37) in the hmit 
Ml = U2 reduce to those of system ( 5.15 ). 

Next we investigate the behavior of the Abehan differentials of the second kind cr^^^ (r, (5, w), fc > 0, 
in the Hmit S —^ 0, when v belongs to one of the gaps ( 5.13 ). We have that JlGf 



crf+^(r, S,v) ~ cr^(r) 



1 



log (5 



-{v) 



(5.38) 



where w!^{r) has been defined in ( 2.14 ) and Q [v) — (u,±/i(w)). From ( 5.38 ) we can obtain the 



expansion for the polynomial P^^^ {r, S,v) = fl{r) 



dr 



namely 



log (5 



(r - v)^j.{r)uj3^{r) 



-{v) 



4ia 



(5.39) 



so that 



Pt\r,6) 



=F log 6 



JQ-(v) ^k\SI JQ-{v) ^Ol?; 



Using (O), (|4.29|-^!30|) and (O^) we have that 



(5.40) 



k=l 



fQ^i-") aft) 
Jq-{v) '^k\'i) 



JQ-(v) '^oUJ 

g+1 

lX\{v, V, u) ^(2/ - l)f^^i_,(u) ;yg^ I ^^l,2. 

JQ-{v) ^oUi 



(5.41) 



1=1 



Adding and subtracting 



29+1 nQ+(v) 



E 



Km 



-(v) i-Ur 



du,qg{u), 



to (5.41) and using ( |2^ ) and ( p^ ) we obtain 



Ui=U2 JQ (v) 



(5.42) 



where $^(r;ti) and W{r) have been defined in ( |5.8| ) and (5.2) respectively. We observe that 



-{v) 
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for all V belonging to the gaps ( 5.13| ) . 

We need to consider also the quantity in the first equation of (5.37) namely 



-1 



U1—U2 



s=o 



Jog S Pl+^ {v + ^/6, 5) log 5 Pl+^ {v ~^/5,5)^ 



5=0 



k=l 



1=1 



' 1 _ P!^l[v-4l.S) \ 

^ 2^5 log 5 \ Pt' (v + VS, S) Pt^ {v-VS,S)) 



<5=0 



9+1 



q,^Av,v,u)y2L (^2V^log4^o'+'(- + V^,<5) Pt\v-^5,5)) 

Using (2.21), ( 5.10| ), (5.40) and (5.31) we simplify the above relation to the form 



5=0 



2^^{v)d.^'^{v;u) + ^gf(«)^(2; - l)ft,(u) I Pf_,iv) P^i^ 



M-)/Q^-7:^f!(e) 



+ 2fiiv)^i'iv)-M^)- 



k=l 1=1 

P^{v,u) 



rQ+iv) s/^N 
jQ-(v) ^ay'i) 



Using repeatedly ( ^.10|) we can write the above relation in the form 



'R9+^{ui,5) R3+^{u2,S) 



logi^nu^-U2) \Pt\ui,S) Pt\^2.5)^ 



UI—U2 



2^l^{v)^B{lKu) + R3{v) - 



Pliv) 



(5.43) 



where $^(r;u) and R^{r) have been defined in (5.S) and ( p^ ) respectively. 
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From ( |5l0| ), and ( p3| ) system can be reduced to the form 



= 2/x2(i;)$f(w; u) + Z3{v) + 



= 



Pf(r) i?f(r)" 



, z = 1,..., 2.9+1 



(5.44) 



where the polynomial Z^{r) has been defined in (p.2|). From proposition 3.3 the last 2^ + 1 equations 



in (5.44) are equivalent to the condition 



Z\r) = 0, 5 > 0. 



Therefore system (5.44) is equivalent for g > to the system 
= $s(w;ti) 

Jq-{v) 



(5.45) 



, i^l,...,2g + l. 



Because of (5.11), when v € (u2j, i*2j-i), 1 < J < .g + 1, it2g+2 = a, we can split the integral of the 



second equation of (5.45) in the form 



Q+Ci") ru2j-i rQ+iv) 

=2 / fii^)<^3ii;u)dC 

J U2j-1 



Therefore ( 5.45 ) is equivalent to ( 5.15| ) for g > 0. 



For (7 = substituting (3.2) in (5.44) it is easy to check that we obtain a system equivalent to (5.15) 



□ 



Proof of Theorem 5.7 



The equations describing the point of gradient catastrophe of the g-j>hase solution can be obtained 
either considering the limit of the {g + 1)— phase solution when three Ricmann invariants coalesce, or 
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supposing that one of the 2g + I distinct Riemann invariants of the g-phase solution has a vertical 



inflection point for t > 0. For proving Theorem 5.7 we follow the latter possibility. 



On the solution of (3^) dxUi{x,t) = {dui{—Xit + wi)) ^ therefore a point of gradient catastrophe 
of the g-phase solution is determined by the system 



du,i-Xit + wi) = 



(5.46) 



Pf(r) Rs{r) 



z = l,...,2g+l 



where 1 <l <2g + 1. We show that system ( 5.46 ) is equivalent to system ( 5.16 ). For the purpose we 
compute explicitly the derivative dui {—)^it + w/). From a generalization of a result in |^ we obtain: 

d Plim) 1 d P|(r) ^^^^^^ 



so that 
dui 



{-\i{u)t + wi{u)) = -6t 



d Pf{r) 



dr P^ir) 



E^ir {U^ - Urn) 



Y\^^=l{ui - Urn) 



llta^tiui - Um) 



-duiqg{u)du,P^{ui) 



^ ^ 9m {^"^y^ 7n — r. 
i^—'f'i ni—n 



9 p9 9 

ra—n ) • 

n—1 ^ m—n 

(5.48) 



In the above relation we need to compute the derivative 

dr^.P^iui) = drP^{r)\r^u,+du,P^{r)\r=u„ 

where Pgir) = + a^r^^^ + • ■ ■ + a^. For computing the derivatives of the normalization constants 



, in P^{r) we need the following proposition. 



Proposition 5.15 [l^J Let uji{r) and uj2{r) two normalized Abelian differentials on Sg. Let 6e ^ = —r= 

Vr 

the local coordinate at infinity and 

k k 

Define the bilinear product 



k>0 
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then 



OUt 



(5.49) 



where ReSr^u \ ^ is the residue of the differential \ ^ evaluated at r = Ui. 
dr dr 

Applying the above proposition to ctq £ind ak, k ~ 0, . . . g—lwc obtain after non trivial simplifications 



/ a? \ 



_d_ 

dui 



a 



9-1 



V "9 / 



/ 1 \ 

ui 

vrj 



( 
1 








g-2 <7-3 



/ fo 

fi fo 

rg-2 Fg-S . . . 

Vrg_i rg_2 • ■ • 





1 



V < J 



(5.50) 



... 

fo 
fi To J 



{2g~i)Pl^^{ui) 
\{2g-l)Pl_,{ui)) 



where the f fe's have been defined in ( 2.19| ). From the above formula we obtain 

du,Pi{ui) ^drP^{r)\r=u,+du,P^{r)\r^u, 



"2"^ °^ ^'-'^2nr^-^,:,,,(..-..),i^, 



£(2n-l)P^i(«0 E<f™_„. 



(5.51) 



Using the relations ( |4.2(]| ), ( [4.24[ ), ( ^.47] ) and ( p.5lD we simplify ( ^.48D to the form 

^2g+l T-r2g+l 



_5_ 

5u/ 



{-\i{u)t+wi{u)) 



2(Po«(^.0)2 °^ 



k^l m^k,l 

WW) 



(5.52) 



+ ^JW^-^9r ( - W E 9™r™-„ - xPo^(r) - 12iPf (r) ) 

^ '/ \n=l m=n / 

where the polynomial Z^{r) has been defined in ( ^.2| ). Applying in the second term of ( 5.52| ) the 
relations (5.1C) we can rewrite (5.52) in the compact form 



— (-AKu)t + z..(.)) = --^ 



WW) 



(5.53) 



where 3>^(r; m) has been defined in (5.8). From proposition 3.3 the last 2g + 1 equations in (5.46) are 
equivalent to the condition 



Z3{r) = 0, g>0. 
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Therefore we can simplify (5.53) to the form 



— {-Xi{u)t + wi{u)) 



711^1 



(5.54) 



when ui > U2 > • • • > U2g+i satisfy the g-phase Whitham equations. As regarding the second 



derivative ■^-^{—Xi{u)t + wi{u)), from ( 5.54 ) we obtain 



—^{-Xi{u)t + 'wi{u)) = — 
duf 



d_ 

dui 



V 



$^(-u;;w) 



-du,'^'{ui-u). (5.55) 



Observing that 

9„,$f(w,;M) =a^$f(r;?r)|,,^„, + a„,$9(r;w)|^^„,, 9^$^ (r; = 2a„,$9(r; 'a)|^^„, , 



we obtain the relation du,^^{ui:,u) = -dr^^ {r;u)\r=ui- Therefore 



-^{~Xi{u)t + wi{u)) = — 



d_ f nm=l,™5^/(w; -Wm) 

dui 



n2g+l( 
[Ui - Ur^ 



From ( ^.54|^.56D and the fact that 



system ( 5.46 ) is equivalent to system ( 5.16 ). 



PoM 



P^M 



(5.56) 



7^ for ui > U2 > ■ ■ ■ > U2g+i, it is clear that 

□ 



We remark that to avoid higher order degeneracies in system (5.16), we impose the condition 

(9,)2$9(r;,?),^„,(,^^,^) ^0. (5.57) 

Indeed we prove that the above condition guarantees that the genus of solution of the Whitham 
equations increases at most by one in the neighborhood of the point of gradient catastrophe. It is 
sufficient to show that the transition to genus g + 2 does not occur. For the purpose let us suppose / 
even and let us consider the function 

^ (9^)^$^ (r;ti)r^„,(^^_4^) = const x $f+^(u/;Mi, . . . , m;, it;, u/, u/, u/, . . . ,^25+1)- (5.58) 



By proposition 5.3 in order to have a genus g + 2 solution in the neighborhood of the point of gradient 
catastrophe, the function 



$^+^(r;ui, . . . + ei,ui + e2,u/ + e3,ui + £4,?/; + £5, ■ ■ .,^23+1) 
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has to change sign in each of the intervals (u/ + £2, + 63) and (u; + £4, + £5) (I even) for arbitrary 



small £1 > £2 > £3 > £4 > £5 > 0. Because of ( |5.58| ) 

^3+'^{r;ui, . . .,ui^i,ui + ei,ui + £2,u; + £3,^; + £4,"; + £5,'«i+i, 



for sufficiently small £1 > £2 > £3 > £4 > £5 > and for r G (u; + £2, + £5). Therefore the transition 
to genus (7 + 2 does not occur. We can exclude transitions from genus g to genus g + n, n > 2, by 
perturbations arguments. Therefore it is legitimate to consider the point of gradient catastrophe that 
solve ( ^.16 ) and satisfies ( 5.57 ) as a point of the boundary between the domains Dg and Dg+i. We 
remark that in the genus 5 = case, the condition (5.58) is not essential as illustrated by theorem 



(3.2) 



Proof of Theorem ( |5.9D 

A phase transition may occur between the zero-phase solution and the two-phase solution. As shown 
on Figure 



5.5 



we can have a double leading edge, a trailing-leading edge and a double trailing edge. 
There are also other types of boundary between the zero-phase solution and the two-phase solution 
which we call "point of gradient catastrophe & leading edge" and "point of gradient catastrophe & 
trailing edge" . Multiple transitions may also occur between the g-phase solution and the {g + 2)-phase 




Double leading edge Leading-trailing edge Double trailing edge 



Figure 5.5: Three kinds of phase transition from the zero-phase solution to the two-phase solution 

solution. In order to determine the systems which describe such phase transitions we consider the 
Riemann surface Sg+2 of genus g + 2 defined by 

jl^ — [r — vi ~ \/Ji)(r ~ vi + V5i)[r — V2 — VS2){r — f 2 + V^2)m^, w G IR , 



2g+l 

n 



fi^ir) = H ~ "j): Ul > U2 > • • • > U2g+1, 



where Vj ^ Ui, j ~ 1,2, i — 1, . . . , 2g + 1. Then we study the hodograph transform for the distinct 
variables vi ± v^i, V2 ± v^2 smd ui, . . . , U2g+i in the independent limits (5i — > and S2 — *■ 0. When 
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vi and V2 belong to the bands (5.12), we are considering the double leading edge. Repeating the 
calculations done for the single leading edge we can determine the equations which describe the phase 
transition for the double leading edge, namely 



= 

i^>9{v2;u) = 



(5.59) 



-I2t 



1,...,25 + 1, 5>0 



From ( ^.14 -5.15) analogous systems can be obtained for the trailing- leading edge and double trail- 
ing edge. 

For studying the phase transitions between the g-phase solution and the (g -I- n)-phase solution, rt > 1, 
having ni leading edges and ^2 trailing edges, rii -t- n2 = n we consider the Riemann surface 



J=l k=l 



25+1 



'Ij), Ui> U2> ■■■> U2g+l. 



Here Vj, j = 1, . . . belongs to the bands (5.12) and y^, fc = 1 . . . ,77.2, belongs to the gaps (5.13) 
< (5j <C 1, j = 1, . . . , ni, and < <C 1, fc = 1 . . . , ^2. 

We study the hodograph transform ( ^.6[ ) for the variables Vj ± j — ^ 
k ~ 1 . . . ,n2 and ui, . . . , U2g+i in the independent limits Sj 0, j = 1 . . . , ni and pk 



■ ,"1, Vk ± y/P^ , 

fc = 1 . . . , 712. 

Repeating the calculations done for proving Theorem and Theorem 5^ it is easy to show that 
system ( ^.23 ) describes the phase transition between the g-phase solution and the (5+7T,)-phase solution 
having ni leading edges, n2 trailing edges and no points of gradient catastrophe, ni + n2 = n > 1. li 
we suppose that the (7-phase solution has also points of gradient catastrophe, 71,1+712+77-3 = 77, 
then combining Theorem 5.4, Theorem 5.(: and Theorem 5.7 we obtain proposition 5.12. □ 



6 Conclusion 



In this work we have constructed, in implicit form, the g-phase solution of the Whitham equations for 
monotone increasing smooth initial data x = f{u)\t=o. The goal is obtained solving the Tsarev system 



(3.4). We have shown that the solution of the Tsarev system which satisfies the natural boundary 
conditions (3.7- p.lO ) is unique. Then we have investigated the conditions for the solvability of the 
hodograph transform (^.6[) . For the purpose we have derived all the equations which describe a phase 
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transition of the solution of the Whitham equations. Studying when phase transitions occur we have 
been able to prove the second main result of this work. Namely we have shown that when the initial 



data satisfies (5J.), the solution of the Whitham equations has genus g < N for all x and t > 0. It is 
still an open problem to effectively determine, on the a; — t > plane, the function < g{x, t) < N 
from the generic initial data x — f{u)\t=o- 
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